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ABSTRACT 

We make first report on inverse chaos synchronization between bi-directionally non-linearly and 
linearly coupled variable multiple time delay Ikeda systems. The results are of certain importance 
in secure chaos-based communication systems. 
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1. INTRODUCTION 

In recent years the presence of chaoic vibrations and its control in nonlinear dynamical systems 
in physics, power electronics, ecology, economics and so on, has been extensively demonstrated 
and is now well established,see,e.g. references in [1-2]. Chaos synchronization [3] is another basic 
feature in nonlinear science and is one of fundamental importance in a variety of complex phys- 
ical, chemical, and biological systems, see e.g. references in [1-2]. Potential application areas of 
chaos synchronization include secure communications, optimization of nonlinear system perfor- 
mance, information processing, and pattern recognition phenomena [1-2]. 

Recently, delay differential equations [4] have attracted much attention in the field of nonlinear 
dynamics. The high complexity of the multiple time-delayed systems can provide a new architec- 
ture for enhancing message security in chaos based encryption systems [5] . In such communication 
systems message decoding would require chaos synchronization between multiple time-delayed 
transmitter and receiver systems [1-3]. Variable multiple time delay systems are further general- 
ization of the fixed time delay systems [6]. Such variations of time delays could be intentionally 
or as a result of fluctuations. In a word, modulated time delay systems could be more realistic 
models of interacting systems. Investigation of synchronization possibilities in such systems are 
of certain importance. 

There are different types of sychronization in interacting chaotic systems [1-2] . Complete, gen- 
eralized, phase, lag and anticipating synchronizations of chaotic vibrations have been described 
theoretically and observed experimentally. Complete synchronization implies coincidence of states 
of interacting systems, y(t) = x(t) [3]; a generalized synchronization is defined as the presence 
of some functional relation between the states of response and drive, i.e. y(t) = F(x(t))\ phase 

1 



synchronization means entrainment of phases of chaotic oscillators, n$ x — m$ y = const, (n and 
m are integers) whereas their amplitudes remain chaotic and uncorrelated; lag synchronization 
appears as a coincidence of shifted-in-time states of two systems, y(t) = x T (t) = x(t — r) with 
positive r; anticipating synchronization also appears as a coincidence of shifted-in-time states 
of two coupled systems, but in this case, in contrast to lag synchronization, the driven system 
anticipates the driver, y(t) = x(t + r) or x = y T ,r > [1-2]. An experimental observation of 
anticipating synchronization in external cavity laser diodes has been reported in [7]. For inverse 
synchronization [8] , a time-delayed chaotic system x is coupled with an another system y in such a 
way that one system's dynamics synchronize to the inverse state of the other system or vice- versa: 
x{t) = -y(t). 

Most chaos based communication techniques use synchronization in unidirectional master-slave 
system. Such a coupling scheme prevents the messages being exchanged between the sender and 
receiver. A two way transmission of signals requires bidirectional coupling. With this in mind this 
paper presents the first report of the inverse chaos synchronization between bi-directionally non- 
linearly and linearly coupled modulated multiple time delayed Ikeda models with two feedbacks. 



II. SYSTEM MODEL 

First consider inverse synchronization between bi-directionally nonlinearly coupled variable time- 
delayed Ikeda systems with two feedbacks, 

dx . . T , . /1N 

— = —ax — mi sin x Tl — m 2 sin x T2 + K y sin y T[i (1) 

(Jbv 

dy 

— = -ay - m 3 sin y T1 - m 4 sin y T2 + K x sin x T3 (2) 
at 

This investigation is of considerable practical importance, as the equations of the class B 
lasers with feedback (typical representatives of class B are solid-state, semiconductor, and low 
pressure C0 2 lasers) can be reduced to an equation of the Ikeda type, see e.g.references in [8]. The 
Ikeda model was introduced to describe the dynamics of an optical bistable resonator, playing 
an important role in electronics and physiological studies and is well-known for delay-induced 
chaotic behavior, see e.g.references in [8]. Physically x is the phase lag of the electric field 
across the resonator; a is the relaxation coefficient for the driving x and driven y dynamical 
variables; r 12 = T"oi,02 + Xi(t)T a i >a 2 sin^i^) are the variable feedback loop delay times;r 3 = 
t~03 + Xi(t)T a 3 sin^t) is the variable time of flight between systems x and y;Toi i o2,o3 are the zero- 
frequency component, r a i ia2i a3 are the amplitude, are the frequency of the modulations; x\{t) 
is the output of system (1) (K y = 0) for constant time delays, i.e. T\ = t 01 ,t 2 = r 02 ; m 12 and 
m 3 4 are the feedback strengths for x and y systems, respectively; K x>y are the coupling strengths 
between the systems. Variable time delays 71,2,3 are chosen in such a way to include both chaotic 
and non-chaotic components. 
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Before considering the case of modulated time delays, we present the conditions for inverse syn- 
chronization for fixed time delays, i.e. 0^2,3 = 0. 

One finds that systems (1) and (2) can be synchronized on the inverse synchronization regime 

x = -y (3) 
as the synchronization error signal A = x — (—y) = x + y for small A under the condition 

m 1 = m 3 , m 2 = m 4 , K x = K y (4) 

obeys the following dynamics 
dA 

= —a A — m 3 A T()1 cosx roi — m 4 A T02 cosx T02 + K y A T03 cosx T03 (5) 

It is obvious that A = is a solution of system (5). It is noted that solution x = —y under the 
conditions (4) is also follows from the symmetry of Eqs.(l) and (2). 

One can study the stability of the synchronization regime x = —y by using the Krasovskii- 
Lyapunov functional approach. According to [4], the sufficient stability condition for the trivial 
solution A = of time-delayed equation ^ = — r(t)A + Si(i)A Tl + s 2 (t)A T2 + s 3 (t)A T3 is: r(t) > 
|si(t)| + |s 2 (t)| + |s 3 (t)|. 

Therefore by using the Razumikhin-Lyapunov functional approach we obtain that the sufficient 
stability condition for the synchronization manifold x = —y can be written as: 

a > |m 3 | + |m 4 ] + \K y \ (6) 

As Eq.(5) is valid for small A stability condition (6) found above, holds locally. Conditions (4) 
are the existence conditions for the synchronization regime (3) between bi-directionally coupled 
Ikeda systems (1) and (2) with multiple delays. 

At first glance, the sufficient stability condition (6) for chaos synchronization is difficult to 
satisfy, as higher values of a could render the dynamics trivial, i.e. non-chaotic. Fortunately 
there is a way out of this impasse. The key point is that the stability conditions derived from the 
Lyapunov-Razumikhin approach is a sufficient one: it assures a high quality synchronization for 
a coupling strength estimated from the stability condition, but does not forbid the possibility of 
synchronization with smaller coupling strengths. Below we demonstrate numerically that one can 
still achieve high quality chaos synchronization in mutually coupled Ikeda systems with double 
time delays without the condition (6) being fulfilled. 

In the case of variable time delays establishing the stability conditions for the synchronization 
is not as straightforward as for the constant time delays. Having in mind that for ui = we 
obtain a case of constant time delays, then as an initial guess one can benefit from the existence 
conditions for the constant time delays case. It is our conjecture that high quality synchronization 
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x = —y could be obtained if the parameters satisfy conditions (4). As evidenced by the numerical 
simulations below, this conjecture is found to be well-based. 
Now consider linearly bidirectionaly coupled Ikeda systems: 

^ = —ax — mi s m x n — m 2 sin x T2 + Ky (7) 

(JjL 

dy 

— = —ay — m 3 sin y T1 — m 4 sin y T2 + Kx (8) 

UjV 

For this kind of coupling the existence conditions for inverse synchronization regime (3) are: 
mi = m-s and m 2 = m 4 ; with sufficient stability condition: a > |m 3 | + |m 4 | + K 

III. NUMERICAL SIMULATIONS 

In numerical simulations to characterize the quality of synchronization we calculate the cross- 
correlation coefficient C [9] 

C ( A q = < (x(t)- < x >)(y(t + At)- <y>)> 

^J< (x(t)- < x >) 2 >< (y(t + At)- < y >) 2 > ' 

where x and y are the outputs of the interacting laser systems; the brackets< . > represent the 
time average; At is a time shift between laser outputs. This coefficient indicates the quality of 
synchronization: C = —1 means perfect inverse synchronization. 

As mentioned above, in chaos based communication schemes synchronization between the 
transmitter and receiver systems are vital for message decoding. With this in mind in the remain- 
der of the paper we focus on the inverse synchronization between bidirectionally non-linearly and 
linearly coupled Ikeda systems with double variable time delays. Figure 1 portrays time series 
of the Ikeda system x (solid line) and system y (dotted line) for inverse chaos synchronization 
x = — y between non-linearly coupled systems, Eqs.(l) and (2) for variable feedback time de- 
lays Ti(t) = 3 + 2xi(t) sin(0.15t), r 2 (t) = 5 + 2x±(t) sin(0.15t) and variable coupling time delay 
Ts(t) = 7 + 2x\(t) sin(0.15t) with parameter values as a = 3, m\ = m 3 = 3.1, m<i = m 4 = 2.5, K x = 
K y = 0.03. xi(t) is the solution of Eq.(l) (K y = 0) for r i = 3 and r 02 = 5. Figure 2 depicts 
synchronization error dynamics A = x + y versus time for parameters as in figure 1. C =0.99 is the 
cross-correlation coefficient between the transmitter x and receiver y system outputs. It is noted 
that the parameters used here in the numerical simulations satisfy the existence condition (4), but 
fail to satisfy the sufficient stability condition (6). Nevertheless high quality synchronization is 
achieved. Indeed, the value of the cross-correlation coefficient C testify to the high quality chaos 
synchronization, which is vital for information processing in chaos-based communication systems. 

Next we present the results of numerical simulations for the linearly coupled Ikeda systems. 
Figure 3 presents time series of the Ikeda system x (solid line) and system y (dotted line) for in- 
verse chaos synchronization x = —y between linearly coupled systems, Eqs.(7) and (8) for variable 

4 



feedback time delays Ti(t) = 1 + 0.5xi(t) sin(O.li), r 2 (t) = 3 + 0.5xi(t) sin(0.1£) with parameter 
values as a = 4, mi = m 3 = 3.2, m 2 = m 4 = 4.5, K = O.Ol. Xi{t) is the solution of Eq.(8) (K = 0) 
for r i = 1 and r 2 = 3 

Figure 4 shows receiver output y versus transmitter output x for parameters as in figure 3. C =1 
is the cross-correlation coefficient between the transmitter x and receiver y system outputs. As 
in the case of non-linear coupling despite the failure to satisfy the sufficient stability condition 
a > \m^\ + |?tt,4| + K, as testified by the value of cross-correlation coefficient (C = 0.99) we were 
able to achieve high quality chaos synchronization. 

IV. CONCLUSIONS 

To summarize we have reported on inverse chaos synchronization in bidirectionally non-linearly 
and linearly coupled variable multiple time delayed Ikeda systems. These results are of certain 
practical importance in secure chaos-based communication systems. 
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Figure Captions 



FIG.l. Numerical simulation of bidirectionally non-linearly coupled variable time delay sys- 
tems, Eqs.(l-2) for a = 3, m\ = m 3 = 3.1, m<i = 771,4 = 2.5, K x — K y — 0.03 and T\{t) = 
3 + 2xi(t)sin(0.15t),r 2 (t) = 5 + 2zi(t) sin(0.15t), r 3 (t) = 7 + 2xi(t) sin(0.15t). x^t) is the solution 
of Eq.(l) (K y = 0) for T01 = 3 and tq2 = 5. Inverse synchronization: Time series of x system 
(solid line) and y system (dotted line). Dimensionless units. 

FIG. 2. Numerical simulation of bidirectionally coupled variable time delay systems, Eqs.(l-2). 
Error dynamics, A = x + y versus time. The parameters are as in figure 1. C is the cross- 
correlation coefficient between the Ikeda systems. Dimensionless units. 

FIG. 3. Numerical simulation of bidirectionally linearly coupled variable time delay systems, 
Eqs.(8-9) for a = 4, mi = m 3 = 3.2, m 2 = m 4 = 4.5, K = 0.01 andn(t) = l+0.5:ri (t) sin(O.lt), r 2 (t) = 
3 + 0.5^1 (t) sin(O.lt). xi(t) is the solution of Eq.(8) (K = 0) for t i = 1 and r 2 = 3. Inverse syn- 
chronization: Time series of x system (solid line) and y system (dotted line). Dimensionless units. 

FIG. 4. Numerical simulation of mutually linearly coupled variable time delay Ikeda systems, 
Eqs. (7-8). Correlation plot between x and y. C is the correlation coefficient between x and y. Di- 
mensionless units. 
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